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i^»» , Recent local tunneling data indicate strong nanoscale inhomogeneity of superconducting gap in 

^vj ' high temperature superconductors. Strong local nanoscale inhomogeneity in the bosonic scatter- 

ing mode has also been observed in the same samples. We argue that these two inhomogeneities 
^y, directly related to each other. To address local boson scattering effects, we develop a local strong 

coupling model of superconducting pairing in a coarse grained superconducting state. Each patch is 
characterized by local coupling to the bosonic mode as well as by local mode energy. We find that 
local gap value on each patch grows with the local strength of electron-boson interaction. At the 
Cm , same time local gap value decreases with the local boson mode energy, an observation consistent 

with the tunneling experiments. We argue that features in the tunneling spectrum due to boson 
(-h'' scattering are consistent with experimentally observed spectra. We also address the ^^O to ^^O 

Q , isotope substitution. Since both coupling constant and boson energy could change upon isotope 

^ ' substitution, we prove that interplay between these two effects can produce results that are very 

I ' different from conventional BCS model. 
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C^ . I. INTRODUCTION 



Electron-boson interaction is at the center of the pairing interaction in conventional superconductors. Pioneering 

work of Eliashberg^ McMillan and Rowell;^, and Scalapino^ conclusively proved that the non-BCS features in the 

tunneling spectra in conventional superconductors are directly related to the electron-boson coupling and ultimately to 

^ ■ the formation of the superconducting state. Experimental evidence of the "strong-coupling" features^ in the tunneling 

spectra was clearly connected to the known phonon spectra in these materials^ On the other hand the number of 

_. 1 ' superconducting materials where electron-boson coupling has been seen was small and initial weak coupling approach 



of Bardeen-Copper-Schrieffer theory was successful in predicting all the measured properties in these materials. The 
few superconductor, like Pb and Sn, with clear electron-phonon features in the spectrum were called "bad actors" 
^\ ' because of deviations from BCS predictions. 

C^ . In this paper we will focus on high temperature superconductors (high-Tc). In high- Tc materials situation is very 

\l ' different. Tunneling spectra in all of these materials clearly deviates from the mean field BCS d-wave tunneling density 
of states (DOS), see e.g. Fig.Q. In this sense all of the high-Tc materials are "bad actors". 

There has been substantial evidence for strong quasiparticle renormalization from tunnelingSiLSiSiifliiiiiSii^ 

and angle-resolved photoemission (ARPES)-^^i-^^i-^^i-^'''i-^^i-^^i^'-'i^^i^^i^'^i^^i^^i^^i^''' experiments. The strong quasi- 

"jrt ' particle renormalization has been suggested as a manifestation of strong electronic coupling to collective 

Cj [ modesj-^^i^^i^°i^-^i^^i^'^i^"^i^^i^^i^^i^^i^^i'^''i^-^i^^i^^i^'^i^^i^^i^'''i'^^ Observed strong electron-electron correlations are clearly im- 

I . portant for mechanism of superconductivity. At the same time the physics of electron-boson coupling in these materials 

'^ ■ is crucial for understanding of the pairing in these materials. We would argue that these features of the tunneling 

G , spectra are due to strong electron-boson coupling. 

>^ ■ Recent scanning tunneling microscope (STM) experiments reveal strong electronic inhomogeneity^2ii2i4iii2ii2ii4 

.. [ More recent experiments with inelastic electron tunneling spectroscopy (lETS STM) have allowed one to directly 

^ ' observe inelastic tunneling features in Bi2Sr2CaCu208+5 (BSCCO) high-Tc materials. These experiments have shown 

k>( [ that bosonic modes that produce strong coupling features in tunneling spectra are also inhomogeneous on the scale 

'V^ . of 20 — bOAm^ Nontrivial correlations between local gaps and local boson mode energies were observed. 

d ' Here we follow the notion that deviations from the BCS tunneling DOS in these materials are caused by the 

strong coupling effects due to electron-boson interaction in these materials. We present a strong coupling theory of 

d-wave superconductor where superconducting state arises as a result of a pairing mediated by bosonic modes, that 

is attractive in d-wave channel. 

In this regard; the treatment is analogous to the conventional Eliashberg-McMillan-Rowell approach. The main 
departure from the conventional approach is that we explicitly allow inhomogeneity in the electron-boson coupling 
strength and the bosonic mode energy. Typical size of inhomogeneity we will assume is on the order of 20 — 50A. We 
will assume that the pairing is local and determined by local values of coupling strength and mode energy at the given 
patch. This approximation allows us to simplify the calculation dramatically. If one takes guidance from the data, it 
is clear that inelastic tunneling features and superconducting gaps are rather local, and there is no "self-averaging" 



seen in the tunneling spectra. In other words, local approximation might be a good starting point for this kind of 
analysis. 

We solve self-consistently the Eliashberg equations on each patch and find local d-wave order parameter. Random 
Gaussian distributions of the local coupling constant and local mode energy are considered. As the result, local 
order parameter (OP) values are random maps that correlate with the input parameters. We find that the local OP 
positively correlate with the coupling constants. One of the most important findings of this local formalism is that 
we find indeed anticorrelation between the typical local boson frequency and the superconducting OP. This negative 
correlation is a direct consequence of the strong local coupling nature of the pairing we assumed. The lETS-STM 
experiment by Lee et alm^ has indeed shown a negative correlation between gaps and lETS mode energy. 

Using our local model we also address isotope effect. The most commonly used isotope substitution for high- 
Tc materials is ^^O by ^^O substitution. Isotope effect was studied in the past by measurements other than STM 
tunneling. Here we point to the important papers in this regard by Shen and Lanzara groups — that argue how 
isotope effect can change both the characteristic bosonic frequency and the coupling constant to the local bosonic 
modepiSiSLSSiSSiSl We find that isotope shift in these two quantities can either work together to mutually enhance the 
superconducting gap, or they can work against each other partially cancelling and therefore making net isotope effect 
small. Hence we conclude that the naive arguments about the "smallness" of isotope effect with substitution of ^^O 
by ^^O are misleading. 

While in this paper we focus on electron-boson interaction, in real systems this interaction contributes to the 
pairing in these materials. We would like to make it clear that we believe pairing in high-Tc materials is a result 
of interplay between strong electron-electron correlations^ and electron-boson interaction. To address the effects of 
spatial inhomogeneity of tunneling lETS spectra we focus only on electron-boson coupling, ignoring electron-electron 
interaction part that will not produce lETS features. It is not clear how high the transition temperature would be 
assuming only electron-boson pairing. We leave this question for a separate investigation. 

The outline of the paper is as follows: In Sec. II, we will introduce local strong coupling formalism, and outline 
the details of the formalism starting with the general inhomogeneous pairing theory in reals space. In Sec. HI, we 
will present results from solving self-consistently the strong coupling equations. In Sec. IV we will discuss the isotope 
effect for inhomogeneous superconductor. We conclude in Sec. V. 

II. LOCAL PAIRING FORMALISM 

We start with the strong coupling analysis in the local pairing limit. Locality here would be understood in the 
sense of a coarse grain approach where we assume that the typical sizes of grains are on the order of the coherence 
length of superconductor ^ ~ 20A. This assumption is consistent with the STM observed granularity in Bi2212 ,Mi^ 

We will present a detailed description of local strong coupling theory with the steps in the logic outlined. Some 
of the points are well known but we keep them in for completeness. The formulas are very similar to the standard 
Eliashberg discussion except we want to stay with the real space description. iSi^ 

Taking the STM data as guidance we can imagine that we take a coarse grained vies of the sample. We are taking 
the given field of view and pixelizing it in a set of boxes with characteristic size of 20 x 20A. Each of the pixels will 
be assigned three variables: order parameter W^(R), bosonic mode energy r2o(R) and local coupling constant ^(R). 
Coupling constant and bosonic mode energy variables are randomly drawn from the given distribution. We calculate 
the local order parameter self-consistently. 

To start we write down the Green's function equations in the Nambu space, assuming no translational invariance. 
We consider the case of both electron-electron and electron-boson interactions being present. The electron-boson 
coupling term in the Hamiltonian assuming no translational invariance is: 



H, 



e — lattice 



= J drdr'dr"cl,Cr'Abl' + W")9{v, r', r") , (I) 



where ^(r, r',r") is the coupling constant and x{y) = b^ + b]. is the boson field. If we assume that only electronic 
density coupled to the lattice degrees of freedom, then coupling constant g(r,r',r") = g(r)5r,r"^r,r"- The bosons that 
couple to the electronic density on the other hand are taken to be on the bonds connecting nearest neighbor sites, 
or alternatively to reside on the dual lattice. We will focus on the local bosonic mode coupling. This seem to be a 
general enough situation that we believe will capture the relevant physics to be addressed. 
Next we consider the second order terms in the effective action that would look like 

j drr'r"drir'/r'i" j drdngir, r', r").9(ri, r'^, r'/)cL(r)crv(T)ct^,^ {riW^a, (ri)xr"(r)M' (^i) ' (2) 



here we assume that phonon propagator B{r — ri,T — ti) = —{Trx{r){T)x{ri){Ti)) is local given the STM data that 
indicate the strong loeal variations of the bosonie excitations as seen in lETS STM^. Thus we assume 

-(T,x(r)x(ri))=B(r,T)(5r,r' (3) 

Then in the self-energy terms electron-boson interaction would produce the term that we would look as 
g{r)g{r"')BQ^{r,r'\r" ,r"'). We again stress here that both coupling constant g and boson energy will be assumed 
inhomogeneous. This situation is qualitatively different from weak coupling approach where only single effective 
coupling constant will be position dependent. 

Equation for Green's functions in Matsubara frequency a;„ = TT(2n + l)kBT becomes: 

W*{Y,Y',lUJn) -Er- S{r,r\lUJn) J<y V 7 , n; v ; yj 

with (^ understood as a convolution in real space, e^ is the kinetic energy operator, that in momentum space 
will have the form £(k) = — i(cos(fca;a) -|- cos{kya)) — pi, fi being the chemical potential. 1 being the indentity 
matrix, and G{r,r' ,iuJn) being a matrix in Nambu space with relevant components G(r,r',ia;„)ii ~ G'(r,r',iw„) and 
G{r,r' ,iujn)i2 = F{r,r',iujn) being the normal and anomalous Green's functions. 

We also explicitly keep the normal self energy in the Gorkov equations: 5(r,r',ia;„) that renormalizes the normal 
selfenergy propagators. 

Normal and anomalous selfenergies are defined selfconsistently through: 

5(r,r',zo;„) = -T^ /'dr"dr"'K.„_^„(r, r'|r", r"')G(r", r'", zc.™) (5) 

W^(r,r',zc.„) = -T^ /'dr"dr'"K.„_^„(r, r'|r",r"')^(r,r',z.;„) (6) 

m '' 

where the effective pairing interaction Vh„(r,r'|r", r'") is given by the combination of both direct and lattice- 
vibration-induced electron-electron interaction: 

y(r,r'|r",r"')n„ = V^-(r, r'|r",r"') - g(r)g(r"')i?o,„(r, r'|r",r"') (7) 

V is defined on Matsubara bosonie frequency 17„, ~ 2TrinkBT, assuming local coupling g. Here V'^'^ is an electron- 
electron interaction written in real space (below we will assume that this term might be inhomogeneous as well). We 
assume V^"^ to be weakly frequency dependent. 

Next we introduce local (on a coarse grained scale) description for the properties of superconductor. For any 
discussion of the local nature of pairing we will need to keep track of the relative coordinate and center of mass 
coordinates. For example, consider the pairing amplitude W{r,r',ujn) that describes the pairing amplitude of two 
particles at sitesr and r'. It is convenient to introduce the center of mass and relative coordinates for the pairing field 
and for the kernel g-B 

R = l/2(r + r'), Ri = l/2(ri + r'^) 

f = r — r', fi := ri — r'j^ (8) 

for simplicity of notation we will drop sign in f hereafter with understanding that capital coordinate label would 
mean center of mass coordinates and small coordinate label would mean the relative coordinates. In thee coordinates 
W{r,r',u!) becomes W(R,r,uj). And similar expressions for S,A etc. The interaction term we will also have a 
factorizable form 

yn„(r,r'|ri,r'i) = V^o„(R,r|Ri,ri) (9) 

In the case of homogeneous pairing W is independent of R. We will focus on the inhomogeneity of W^(R, r, cu) as 
a function of R. 

We introduce the basis functions for d-wave channel and ignore any other pairing channels. This is not a principal 
assumption but a useful one that allows us to greatly simplify equations. 



We have therefore 

W{R,r,u) = W{R,uj)r]{r) 
ya„(R,r|Ri,ri) = l/a„(R,Ri)??(r)r/(ri) (10) 

and similar for S{r,r',ujn) etc. 

Here ?7(r) is a real space representation of the basis function that has a d-wave character. The simplest way to 
present is is to take a function that is nonzero at nearest neighbors of siter that has a d-wave signature: ry(r) ~ 
^s=±x ±y{~^y^r.r+s on the lattice. In the continuum we would have to deal with the gradient operator: ?7(r) ^ 
(9^ — dy). In the momentum space it will be a simple ?7(k) = cos(fc2:a) — cos(fcj^a). 

It is also convenient to introduce mixed representation where we use Fourier transform for the relative coordinate. 
Then 

WK,R,k) = M^K,R)77(k) 
S'(a'„,R, k) = S'(a;„,R)lk 
Fj2„.(R,k|Ri,ki) = Fn„(R,Ri)r7(k)r?(ki) (11) 

We consider d-wave channel only in assuming a simple factorizablc approximation. This is definitely an oversim- 
plification since the inhomomgeneous system would admit the mixture between components. One can always include 
the mixing in a more detailed approach. In practice this mixture could be modest. For us the main focus here would 
be on the real space modulations of the gap function II^(R, Wn), electron-lattice and electron-electron interactions 
Vh„(R, R')- We proceed with this simplifying assumption that would make our discussion more transparent. 

We can rewrite the selfenergies 5", W Ea. H5lt)|l upon projection on the d-wave channel as: 

W{R,iUr,) = -TY, /"dR'K.„_^„(R,R')(^(a.™,R',k'))d (12) 

and similarly for S: 

5(R,c^„) = -tY, /dR'K.„-c.„(R,R')(G(cc',„,R',k')) (13) 

with 

(F(c^,„, R', k'))d = j dk'77(k')F(c^,„, R', k') (14) 

{S{uj^, R', k')) = j dk'5(a;,„, R', k') (15) 

We focus on the gap equation hereafter. We find that normal self-energy 5'(a;„,R, k), at most leading to mass 
renormalizations on the scale unity. The fermi surface average correction due to normal self energy corrections is 
small and hence we ignore it. This allows us to keep only d-wave projected part of the interaction in Ea. (|ll|l . 

From the solutions of the Green's function we have self consistently defined F, G: 

pf R , ^ ^^„(R)^(k) ,^.. 

^ "' ' ^ ~ {lU:n ~ 5(C.„,R))2 - £(k) - M/2(^„^R)^2(k) ^'^^ 

r(,, -R W^ - iuj„~ S{uJn,'R) f„. 

U^a;„,±t,Kj (,^„_^(^„,R))2_,(k)_VF2(a;„,R),72(k) ^" > 

Gap in the quasiparticle spectrum is determined as 

.. „. Wiu;„,R) 
A(u;„,R) = - — — — (18) 

1 - b[UJn,R)/ii^n 



These equations arc written in general form. We take S = below. 

Equations Eq. H12I13I1()I17I18|) are the main result of this seetion. These equations are quite general and describe 
the inhomogencous superconducting state in the presence of inhomogeneous pairing interaction. 

These equations are similar to the Eliashberg equations considered for a homogeneous superconductor. Here we 
focus on the spatial dependence of the superconducting properties like gap in the spectrum and pairing interaction 
l/a„(R,R'). 



A. Local Approximation 



We can make further progress if we will make some additional assumptions. We will assume that the kernel in 
Ea ll2ll3|) local. Again, this locality should be understood in a coarse graining sense with typical length scale for 
coarse graining to be on the order of superconducting coherence length ^ ~ 20A. This length scale is compatible with 
the observed inhomogeneities in the tunneling gap and bosonic frequency, as imaged with STM4^ 

Local on-site pairing kernel would be incompatible with the d-wave character of the pairing we assumed here. For 
a moment we will focus on the electron-lattice part of the kernel. It contains an effective coupling g and boson 
propagator _B as a single combination we label g^B. We will assume local approximation for both couling constant g 
and bosonic frequency. Thus: 



Fo„(R,R') = y^^(R)<5RR. -5(R)5(R')Bn„(R)5RR' 



(19) 



Where, following standard Eliashberg approacbiiSi^, we will assume that the bosonic spectral density (on the real 
frequency axis) would have a local character: 



/m(.g2s(R,R',n)) = ng^{Ti)[S{n - f2o(R)) - SiQ + no(R))]5RR' 



(20) 




(W(/?,,,«),Q„(7?,,),g(/?,,)) 



FIG. 1: Illustration of our coarse graining approach to the local strong coupling solution of Eliashberg equations. We pixelize 
a finite system with the size of the patch to be on the order of ^. Then we solve the selfconsistent strong pairing equations on 
each patch. The outcome of the solution is the set of local value for d-wave order parameter, that depends on the local value of 
boson frequency and local coupling constant. We present these values as a pixel dependent vector with these three components 

{W{Ri,Lo),no{Ri),g(Ri)). 



With the recent STM experiments we now have an independent experimental measure of the local bosonic energy 
^10(1^) as a function of position and doping^. Typically, these energies are randomly distributed in a sample with 
characteristic variations on the length scale on the order of 20 — 50A and thus are consistent with our assumption of 
locality. The sample averaged frequency Qq = (r2o(R)} is essentially doping independent and is about 52 mcV, while 
distribution ranges between 40-70 meV for the observed bosonic modes in STM experiments, see Fig([3l- Wc do not 
have a similar experimental information on coupling constant. 

In practice, of course, only the total kernel V will enter into the self-consistency equations and we would not be 
able to differentiate the effects of inhomogcneity in the W due to electron-electron versus electron-lattice interaction 
inhomomgeneity. There is one important distinction however between electron-electron vs electron-lattice coupling. 
Electron-electron interaction being essentially frequency independent can not produce features outside the coherence 



peaks. Electron-lattice coupling on the other hand will produce the features in local tunneling characteristics at 
u;(R) = A(R) + r2o(R). Hence the local tunneling characteristics would allow us to extract the inhomogeneous 
values of the bosonic modes, at least in principle. In practice one would have to deal with rather large signal- noise 
uncertainties but the local bosonic energy extraction from the dl/dV data can be done^. 
The local version of equations Ea H12ll3() now would take a form: 

VF(R,w„) = -T^K.„_^„(R)(i^(a.„„R,k'))d (21) 

777. 

and similarly for S: 

S{R,LOn) = -r^K,„_^,„(R)(GK„R,k')) (22) 

rn 

and one can recognize standard Eliashberg equations that are now written locally, patch by patch. Hereafter we 
will approximate ^(R, w„) = (1 — Z^^^)iuJn, i.e. local Z factor normalization only. In practice we know that effective 
mass renormalization in high-Tc materials is not very large and at most Z ^ 2, hence the effects of the normal self 
energy on quasipatrticle dispersion would be minor. It is the pairing interaction contribution from electron-boson 
interaction that we will be paying attention to. These equations in the homogeneous case were well analyzed^'^ 

B. Weak Coupling Approximation 

Now we will consider the weak coupling limit of these equations, namely the case when the pairing amplitude W and 
normal self-energy corrections S are small compared to the typical bosonic energy, f2o = (^(R-))- The most natural 
region to make the weak coupling approximation would be in the overdoped regime where both superconducting 
gap A and Tc decrease with increased doping. There are no competing orders in the overdoped regime that would 
make the analysis more complicated. Competition with the other orders, like charge ordered state and pseudogap 
would make our analysis in terms of a single superconducting gap inaccurate. Therefore the analysis presented below 
assumes that we are dealing with optimally doped to overdoped samples. 

Summation over Matsubara frequency in Eg. H2 1122(1 is treated in a standard way by using Tj^nfi^n) = 
<f dz ta,nh{^) f (z) , P = l/fc^T. This integral is reduced to the integral over spectral density for the effective in- 
teraction that we will assume to be local: 



/oc 
dx 
-oo 



ImV{x) 



ImigliR, R')Bo(R, R')) = g^^M^ ~ l^o(R)) - S{n + r!o(R))<5R,R' 

imV''%R,R',n) = f'="(R)5r3', \n\<nc 

ImVn(R) = ImV^R) - Img'^BniR) (23) 

where we have parametrized the position dependent electron-lattice interaction by local coupling constant g^ and 
local boson frequency f2o(R). Electron-electron interaction is assumed to be frequency independent up to cut-off 
frequency f2c 3> f2o. 

Analysis we will implement here is essentially the same as the one used in standard Eliashberg approach. We use 
the d-wave projected propagator F and integrate over the quasiparticle energies £(k) to find from Ea (|21l22|l : 

WiR,n) = -A^o/.xr.n-tanh(^)^-(^-^^);i^e(( , ^^^'^'^ ).) 

2g2(R)(»o(R)+nO ^ WjRM') . .„,. 

''° (f7o(R) + nr (ny + ^s''''^^ ^{Zn^n'f - w^{r, k, n') ^"^^ ^ ^ 

In the weak coupling limit for small coupling g'^B(ui,R) we can develop a local BCS pairing approximation for this 
local strongly coupled superconductor by approximating W{uj,R) ~ W{uj = O,R)0(r2o(R) — 'j-'),w > and similar 
step like cut off at negative to. 



We will focus on the low energy part of W(R, ^2) that will contain real part of the gap only. Then a) the integral 
over X is trivial since we assumed spectral function ImV to be a delta function in frequency. We assume T ^ and 
hence tanh(^) ~ 1. b) we assume that electron-electron part will produce normalizations on the normal channel and 
also will produce a d-wave pairing. Then we take a low frequency limit of this equation fi -^ and limit the integral 
over ri' over the range f2' < f^o since the kernel is attractive only in this range. The weak coupling limit therefore 
would read as: 



W(R,fi = 0) = / 

JO 



"«(^) 2.g2(R)(f7o(R) + f7') 



cLVl'Nq 



(r2o(R-) + ^')^ -{^Y + i5 



and ultimately we obtain the local version of BCS equation for W^(R) = W^(R, fi = 0): 

/•f2o(R) W(R) 

with 

J ^ 

is a weak coupling limit coupling constant in Eliashberg theory. 

Using the local approximation for the spectral density Ea. (|23|) . we have 

Here we explicitly assume that .g(R) and il(R) are independent random distributions. This assumption is natural 
if we allow these two quantities to be set by local environment in the crystal and we do not assume here that coupling 
constant .g(R) ~ (Af51(R))^/^ as would be the case for quantized extended collective modes. We thus arrive at the 
local gap equation for A(R) = VK(R)/Z(R): 

A(R) = 17o(R) exp(- ^ ) (29) 

7Vo.9e//(R) 

Equation Ea. (|29(l is applicable in the weak coupling limit and therefore can be viewed only as a qualitative result. 
For any distribution of the coupling constants for small enough average value < g^B > there will be regions where this 
coupling constant is not smaller than 1 and hence weak coupling analysis would fail in those regions. Nevertheless it 
is useful in that it allows us to analyze the results of numerical calculations, see below and compare numerical results 
with the locally observed quantities. With this caveat in mind we will consider the implications of the Ea. (|29|l for our 
analysis of the local pairing. 

We find immediately three important consequences of the local pairing approximation: 1) Effective coupling constant 
is a function of the local boson mode energy. Effective coupling constant geff (R) is inversely proportional to the mode 
energy f2o(R)- This result implies that there is a direct negative correlation between local gap and local bosonic mode 
energy. Similar direct negative correlation is observed in the STM experiments on inelastic tunneling spectroscopy^^ 
2) Another implication of this result is that the isotope efffect will affect both the prefactor and coupling constant 
in Ea. (|29|l . This point will be discussed more below. 3)Finally, from this simple equation we can find local effective 
coupling constant Nogeff(R) as: 

iVoge//(R) = - i(^) (30) 

This equation contains two experimentally observable quantities: A(R) and f7o(R)- We therefore can build the 
real space map of effective coupling constant. 



III. NUMERICAL RESULTS AND DISCUSSIONS 



Here we will discuss the numerical solutions of local Eliashberg equation. From Eqs.(21) and (22) together with 
Eqs.(16) and (17), it can be written explicitly as: 



NlP ^-^ ^ n(R, q; iw„ - ii2,„) 

T^(R,io.„) = -^^^^£££i!^l!Iil^W^(R,q;^c.„-^a.), (31a) 

^lP ^-^ tr" n(R, q; iw„ - iVl^n) 



q iilri 



where 



with 



Me//(R,k;ia„) = [Fee(R) -5'(R)S(R,i^m)]^k (32) 

5(R; i^^) = -^^-— ^^-^ - ———^ , (33) 

and 

n(R, k; lujn) = {N„^(R, q; lUJnf - Kk + W^'(R, k; iw™)]'}^/' . (34) 

Note that the strong correlation between electrons themselves can give rise to an effective pairing. In this case, Vee 
is negative. In principle, we cannot exclude the possibility that Vee also becomes inhomogeneous. Here we focus on 
the effect from the electron coupled to local modes and will ignore the contribution from Wee- We mention equation 
for normal self-energy correction for completeness. As was pointed out, we will ignore S. 

We adopt a six-parameter fit to the band structure used previously for optimally doped Bi-2212 systemSf^ having 
the form 

^k — — 2ii (cos kx + cos ky) — 4^2 cos kx cos ky 
—2t^{cos2kx + cos2ky) 
— 4i4(cos 2kx cos ky + cos fc^ cos 2ky) 
— 4i5 cos 2kx cos 2ky — jjl , (35) 

where ti = 1, ia = -0.2749, h = 0.0872, t^ = 0.0938, t^ = -0.0857, and /x = -0.8772. Unless specified explicitly, the 
energy is measured in units of ti hereafter. 

We use the method of Vidberg and Serene to first solve the above coupled equations in the Matsubrara frequency 
space. On Matsubara axis the quantities 5* and W are real. We then do the analytical continuation with the Pade 
approximation to covert them to the axis of real frequency, on which they have a real and imaginary part. Partly 
motivated by the ARPES experimentsiSiiSiSSiSiiSSiSiSliSiSSiSi, we take the averaged frequency of the local bosonic 
modes to be (rig) = 0.3 and the bare electron-bosonic mode coupling constant (go) = 0.5. The temperature is chosen 
at T = 0.01. To be illustrative, we first show the calculation at a specific coarse-grained spatial point. Our calculations 
show that S is negligible and we will ingore it hereafter. 

A. Relation between features of the gap and coupling constant and bosonic energy 

In Fig.(|21l we can see that at low energies the real part is constant and imaginary part rises only after energy 
exceeds the boson energy VLq = 0.3. The first peak of ReW curve and shoulder on ImW curve precisely correspond to 
the boson energy. Features at higher energies correspond to the multiboson processes and are at multiples of JIq. In 
Fig-® we observe the evolution of the inelastic features as a function of coupling constant. Gap function is changing 
substantially even for fixed boson energy. Still there is always a feature at boson energy for all coupling constants. 
For larger energies and larger coupling constants peaks in ReW are getting broader. In Fig.Q the features in the the 
real part of gap function for different values of boson energy are shown. Again the first shoulder albeit of different 
intensity can be seen at energy that exactly corresponds to the boson energy. Feature at E = 0.2 for blue curve is 
very broad, the one at E' = 0.3, 0.6 can be seen for boson energy fio = 0.3 and finally for the green curve one can see 
feature aX E = 0.4. Features at substantially higher energies are likely to be numerical artifacts of our use of Pade 
approximations in analytic continuation. 
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FIG. 2: The real (red line) and imaginary (blue line) parts of the complex d-wave order parameter as a function of energy. 
Here the bosonic mode frequency Slo ~ 0.3 and the electron-mode coupling constant go = 0.5. 
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FIG. 3: The real part of the complex d-wave order parameter as a function of energy for various values of the electron-mode 
coupling constant go — 0.4 (blue line), 0.5 (red line), 0.6 (green line). The bosonic mode frequency is fixed at flo = 0.3. 

B. Correlation between inhomogeneous coupling constant and gap; anticorrelation between bosonic mode 

energy and gap 

Next we consider the effect of electronic inlioniogeneity. For this purpose, we consider two cases, a) coupling 
constant has a spatial distribution and b) local bosonic mode has a spatial distribution. Both of these parameters 
can be position dependent at the same time as we suspect is the case in real systems. Here we want to differentiate 
between the effects coming from coupling constant and effects coming from frequency variations. We assume that 
both distributions are gaussians. 



P{x) oc exp 



ix-{x)r 



2a2 



(36) 



where x represents go or Qq. Throughout the work, we take a = 0.3. 

In FiglO we observe the direct correlation between the strength of the coupling constant and ReW{E ~ 0). This 
direct correlation is natural and to be expected. In Fig. © we find an anticorrelation between the boson mode energy 
and gap energy. The nature of this anticorrelation follows from our assumption on the boson spectral function that 
is peaked at one energy Oq- Indeed from the structure of the pairing kernel one can see that larger boson energy at 
fixed g would lead to lower effective coupling constant, see Ea. ()27|l . Thus we conclude that the anticorrelation is not 
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FIG. 4: The real part of the complex d-wave order parameter as a function of energy for various values of the mode frequency 
fio = 0.2 (blue line), 0.3 (red line), 0.4 (green line). The bare electron-mode coupling constant is fixed at 30 = 0.5. 
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FIG. 5: The spatial distributions of the bare coupling constant (left panel) and the resultant d-wave order parameter at zero 
energy. The bosonic mode frequency is fixed at Qo = 0.3. 

a consequence of the weak coupling analysis but is present in the full numerical solution of the self-consistent gap 
equations. This anticorrclation is directly observed in the lETS STM experiments^^ 



C. LDOS map and the d-wave order parameter map 

We have also calculated the local density of states (LDOS). 



piE) 



Nltt 



^|Img(R,k;£; + iO+)| , 



(37) 



which correspond to the local differential tunneling conductance as measured by the STM experiments. Fig.lO shows 
the local density of states as a function of energy at a selected spatial point corresponding to Fig.l. For comparison 
we also show typical experimental data for STM tunneling density of states and its derivative, FigQ 

We also calculated spatial image of the LDOS at the energy Ag. In the strong-coupling theory, the d-wave order 
parameter is energy dependent. To demonstrate the distribution of gap inhomogeneity from the LDOS distribution 
at the gap edge Ag, we should use the d-wave order parameter at fixed energy, say E =< A >. FigJHJ shows spatial 
distribution of the LDOS for the case of an inhomogeneous distribution of local bosonic modes. 
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FIG. 6: The spatial distributions of the mode frequency (left panel) and the resultant d-wave order parameter at zero energy. 
The bare electron-mode coupling constant is fixed at go = 0.5. Anticorrelation between the mode frequency and gap magnitude 
is clearly seen. Large mode frequency regions (red spots in the left panel) correlate with the small gap regions (dark spots in 
the right panel). Anticorrelation for arbitrary coupling constant is thus verified numerically. In the weak coupling limit the 
anticorrelation follows from Eq. 1271281 . 
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FIG. 7: The density of states as a function of energy. Here the bosonic mode frequency Hb ~ 0.3 and the electron-mode 
coupling constant go = 0.5. Also shown is an experimentally measured dl/dV and cfl/dV^ spectra at some typical point are 
shown. The overall similarity of the spectrum with the one we have calculated is suggestive that the features in the tunneling 
seen at the large bias are indeed consistent with the strong coupling features due to electron-boson interactions. 



IV. ISOTOPE EFFECT IN INHOMOGENEOUS SUPERCONDUCTOR 



One of the most powerful tools to investigate the role of phonons in the pairing is to study isotope effect. In the case 
of conventional supercoductors it was found that isotope subsitution of the lattice atoms would change the transiton 
temperature thus directly indicating that lattice is involved in pairing. By changing the mass of the ions from M to 
M' critical temperature T^ -- Af"" would change by T^/T^ = {M/M')". 

First consider conventional homogeneous superconductors. In conventional superconductors often a ~ 1/2 is 
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FIG. 8; The spatial distributions of the real part of the d-wave order parameter ai E = Ag (left panel) and the local density 
of states at the same energy. The bare electron-mode coupling constant is fixed at go = 0.5. 

observed. This result follows the observation that the effective coupling constant in standard BCS formalism is 
independent of mass M^ Simple arguments show that effective coupling constant in homogeneous case geff = ,j 2 ; 

where C is a constant, M is the ion mass and cj^^ is the average phonon frequency squared. Since uja-v ^ M^^'"^ for 
phonon spectrum regardless of the detailed shape, one finds that the coupling constant in this case is independent of 
mass M . Therefore the only effect of the isotope substitution is on the change of the phonon spectrum and the cutoff 
frequency that is in the prefactor in BCS equation for 



Tc = VLq exp(- 



1 



Nogeff 



(38) 



We thus find the conventional exponent is a = 1/2 that is set solely by the prefactor fio within the BCS theory. 
Situation we consider is very different. As we pointed out, gi,ff{Ti) is made from two random independent quantities, 
f/(R) and rioCR-)- This will lead to a very different isotope effect in this random superconductor. 

Standard isotope substitution experiment in cuprates is replacement of ^^O by ^^O. The changes in Tc produced 
by this substitution are small, a is small and depend on doping levels of the samples. In the optimally doped samples 
a is essentially zero. 

Situation for inhomogeneous superconductors is qualitatively different from conventional BCS case. Inhomogeneous 
superconductor is characterized by two rather than one parameter that enters into the gap equation Ea. H2U|) : one is 
coupling constant 5(R) and another one is a local boson frequency rio(R')- In principle, both random variable will 
change upon isotope substitution. On general grounds, isotopic substitution would change the local environement as 
it affects both in-plane and out of-plane oxygen atoms. Hence, we argue that both the coupling constant and phonon 
frequency are affected by isotope substitution. It would mean therefore that both prefactor and coupling constant in 
the exponent are changed va Ea. H29|) upon ^^O to ^^O substitution. In addition, for the inhomogeneous superconductor 
one has to differentiate between the isotope shift of critical temperature of a sample Tc and the isotope shift of the gap 
A(R). Here we do not address the net shift of Tc as it is determined by phase fluctuations at higher temperatures. 
Instead we focus on shift of local gap A(R) or average gap < A(R) >. 



A(R) = fJo(R)exp(- 



1 



A^0.9e//(R)' 



(39) 



Again for simplicity we will take a weak coupling limit of inhomogeneous Eliashberg equations. Let us take an 
'average" of the equation as an approximate way to discuss the average shifts in rio(R) and in A^o<?e//(R), Ea. lpli^ 



(A(R)) = (fio(R))cxp(- 



(iVo5e//(R)) 



) 



(40) 



It is known^ that the average frequency fio — (^o(R)) shifts from 52 meV to 48 meV upon ^^O by ^^O substitution. 
On the other hand this shift in Q.^ can be offset by a shift in average je// thus producing an zero isotope effect that 
is very different from BCS exponent aBcs — V^- To illustrate how one gets near zero isotope shift, we take that 
effective coupling constant shifts by 
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1 (^"^ 



Nogeff Nog', 



(41) 



eff 



This shift of the effective couphng constant could offset the shift of the prefactor in Eq. H4U|) . The net isotope effect 
wiU be determined by the combined isotope shift of the boson mode energy and effective couphng constant. They can 
mutuahy cancel each other making net effect small, as we suspect is the case near optimal doping. Alternatively, both 
effects can add up to produce large isotope shift. Thus by knowing only the shift of boson energy is not sufficient to 
address the net isotope effect of the gap. 

To illustrate this point we have preformed "numerical isotope shift" experiment within our model, see Figl^. 
In order to model the effect of oxygen subtitution we change random distribution P{^o) of the boson energy Hq. 
We assumed that changing ^^O to ^^O would shift gaussian distribution with the mean values of boson energy 
(51o)oi6 = ^ = 0.3 to (Oo)oi8 = ^(1 ~ 6%) ( upper panel). At the same time coupling constant g can also change 
upon isotope substitution. We consider the shift in the coupling constant g which here is taken to be constant at the 



same time (lower panel). For °0 we take g 



go 



0.5 and for ^°0 we use .go(l — 4%). We find that negative shift 



of boson energy by 6% can be offset by the shift of the coupling constat go and the net effect would be negative shift 
of the gap (W). To address the net isotope effect it is necessary to measure independently both boson energy and 
coupling constat for ^^O and ^^O. At moment there is no independent experimental measurement of the coupling 
constant we are aware of. 
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FIG. 9: Upper panel: histograms of (the boson energy Qq (upper panel) for O^^ (red) and O^* (blue) are shown assuming 6% 
isotope shift. Lower panel: histograms of the calculated gap distribution W for these two boson mode distribution calculated 
with simultaneous coupling constant shift. We assumed that coupling relative constant shift is negative 4%. This turns to be 
sufficient to produce net negative gap shift. 

For completeness we also present the correlation between the gap and the mode energy Oq, Fig ljlOl) . The negative 
correlation between the gap and mode energy is clearly seen for both O^^ and O^^ isotopes. 



CONCLUSION 



In this paper we use strong coupling model for boson mediated d-wave pairing for inhomogeneous superconductor. 
To model inhomogeneous superconductor we consider a coarse grain model with the typical size of the patch on 
the order of superconducting coherence length ^ '^ 20 — 50A. We use patch dependent pairing interaction due to 
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FIG. 10: From the real space maps for mode energy ilo and gap as a function of position we can extract the statistical cross- 
correlation of these two quantities. We find overall negative concave-shaped A-f^o cross correlation for the cases of O^^ and 
O^^ . Isotope shift of mode energy to lower energies is clearly seen on this plot. The average gap value is affected by the shift 
in the coupling constant and in mode energy. 



disordered pairing boson with patch-dependent coupling constant g(R) and boson energy rio(R-)- This local pairing 
produces local superconducting gap as a selfconsistent solution of Ehashberg equations patch by patch. 

We argue that any inhomogeneous theory of strong coupling theory of pairing has to involve at least two independent 
quantities that characterize electron-boson interaction: coupling constant and boson energy. Reduction of the pairing 
theory to an inhomogeneous BCS pairing model does not allow one to distinguish a relative role of the coupling 
constant g(R) vs r2o(R). In BCS-like analysis one deals with the single effective coupling geff{IV) that is combination 
of g(R) and f^o(R), Eq.IlHll- 

We use the local coupling constant and boson energy drawn randomly with the gaussian distributions and calculate 
local gap function W{lli, R) — ReW{u!, IV) + iImW{u!, R). This gap map is then used to compare correlations between 
VF(R), g(R) and fio(R)- We also calculate the local tunneling density of states that is consistent with the observed 
by STM density of states. 

Numerically we find a direct positive correlation between the gap map M^(R) and boson coupling constant map. 
This is not surprising: the larger the coupling the larger is gap function. We also find an interesting and surprising 
result that there is a negative correlation between gap function scale and boson mode energy Og. We give a simple 
interpretation of this negative correlation in the case of the weak coupling analysis. We find that effective BCS 
coupling constant is inversely proportional to the boson energy, Ea. (|28|l . Because effective coupling constant is in the 
exponent of a gap solution Ea. (|29|l . its dependence on r2o(R) is more important than the dependence on frequency in 
the prefactor. Larger frequency boson is less effective in inducing pairing for fixed ^(R) in our model. This result is 
consistent with the experimental observation of the JETS signal, where anticorrelation between gap and boson mode 
energy was observed]^ 

Exact nature of the boson is not important for our analysis except when we discuss isotope effect. Then we expicitly 
assume that boson is a lattice mode and its energy has on average an isotope shift consistent with the isotope shift 
for phonons due to ^^O to ^^O substitution. 

We also consider isotope effect. Wc find that in order to correctly address the full isotope effect one again would 
need to have assess changes in the coupling constant .g(R) and rio(R) as a result of ^^O to ^^O isotope substitution. 
Both shift in the boson frequency and coupling constant contribute to the net isotope effect and wc find that small, 
on the order of ~ 4%, changes of the coupling constant can completely offset the isotope shift of the gap function 
W^(R) caused by the standard isotope effect for the lattice mode r2o(R)- 
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We therefore find that in order to understand the isotope effect in d-wave superconductors we would need inde- 
pendent measurements of the couphng constant map and boson mode energy map for different isotopes. The whole 
notion that a single isotope exponent can characterise the spatially modulated superconducting state as is the case of 
high-Tc materials seem to be too simplistic to address the real situation. It appears one can not make an evaluation 
on the importance of the lattice effects in high-Tc superconductors based on a shift of a critical temperature without 
addressing the changes of the gap, boson modes and coupling constant. 

To address the effects of spatial inhomogcneity of tunneling lETS spectra we focus on electron-boson coupling, 
ignoring electron-electron interaction part that will not produce lETS features. Electron-boson interaction is only 
one contribution to the pairing interactions. Pairing in high-Tc materials is likely a result of interplay between strong 
electron-electron correlations and electron-boson interaction. The realistic magnitude of pairing interaction and how 
large the transition temperature would be by assuming only electron-boson pairing is an interesting question. We 
leave this question for a separate investigation. 
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